Two-Dimensional Crystals as a Buffer Layer for High Work Function Applications: The Case of Monolayer MoO3

We propose that the crystallinity of two-dimensional (2D) materials is a crucial factor for achieving highly effective work function (WF) modification. A crystalline 2D MoO3 monolayer enhances substrate WF up to 6.4 eV for thicknesses as low as 0.7 nm. Such a high WF makes 2D MoO3 a great candidate for tuning properties of anode materials and for the future design of organic electronic devices, where accurate evaluation of the WF is crucial. We provide a detailed investigation of WF of 2D α-MoO3 directly grown on highly ordered pyrolytic graphite, by means of Kelvin probe force microscopy (KPFM) and ultraviolet photoemission spectroscopy (UPS). This study underlines the importance of a controlled environment and the resulting crystallinity to achieve high WF in MoO3. UPS is proved to be suitable for determining higher WF attributed to 2D islands on a substrate with lower WF, yet only in particular cases of sufficient coverage. KPFM remains a method of choice for nanoscale investigations, especially when conducted under ultrahigh vacuum conditions. Our experimental results are supported by density functional theory calculations of electrostatic potential, which indicate that oxygen vacancies result in anisotropy of WF at the sides of the MoO3 monolayer. These novel insights into the electronic properties of 2D-MoO3 are promising for the design of electronic devices with high WF monolayer films, preserving the transparency and flexibility of the systems.


■ INTRODUCTION
Work function (WF, Φ), defined as the minimum energy required for an electron to escape from a material surface, is a critical surface characteristic as it determines the interface properties in electronic applications. The introduction of high WF transition metal oxides (TMOs) as buffer layers for anodes in organic electronics has proven to increase efficiency of devices due to securing the energy level alignment. 1,2 However, further thickness reduction of conventional TMO layers is a difficult task because it has been reported to impact the electronic properties including their WF. 2−6 Thus, the pursuit of novel designs and alternatives to conventional buffer layer materials is the subject of intense research. In this work, we propose ultrathin (<1 nm) two-dimensional (2D) crystalline films for fine-tuning of the substrate WF, therefore downscaling the thickness of currently used buffer layers. We show that the balance between maximal efficiency and minimal thickness with transparency can be found using an ordered crystallographic structure.
We focus on α-molybdenum(VI) oxide (α-MoO 3 ), a van der Waals (vdW), wide-band gap n-type semiconductor with a high WF of 6.9 eV in its bulk phase. 2,6−8 This value is notably higher than most WFs of typical layered materials such as highly ordered pyrolytic graphite (HOPG) and graphene (4.4−4.5 eV) 2,5 or MoS 2 (4.3−5.5 eV). 9,10 In α-MoO 3 , valence and conduction bands are ∼3 eV apart and the material can be easily doped through oxygen vacancies forming MoO 3−x 7,11 −14 with additional bands below the conduction band minimum and a local decrease in the band gap. 15 It has been theoretically predicted that α-MoO 3 retains its bulk properties in the 2D limit, and thus, the electronic structure of 2D α-MoO 3 is not significantly altered from the bulk, 12,16−19 including its high WF. This is different to, for example, graphene and monolayer MoS 2 , which may present different electronic properties depending on a number of layers. 10, 20,21 Concerning MoO 3 and investigations performed by ultraviolet photoemission spectroscopy (UPS), there are many relevant reports on both the thickness dependence and the influence of the substrate on its WF. 2−4,6,22−24 Much effort has been put into enhancing MoO 3 WF by increasing film thickness or restoring its high value after air exposure. [2][3][4]7,8 The WF of MoO 3−x films was reported to gradually increase with thickness, for example, at 1.4 nm, it reaches 5.90 3 −6.40 eV, 24 while final saturation (≤6.80 eV), that is, the highest observed value, on graphene occurs for thicker films of around 4.0−5.0 nm. 2−4 Although successful, these studies did not investigate the electronic structure dependence on the number of layers of α-MoO 3 and did not focus on potential improvement by full crystalline coverage essential to accurately describe MoO 3 monolayer properties. The high WF of 2D-MoO 3 has not been fully experimentally revealed in part due to the uncontrolled environment and underestimation of results from conventional analytical tools.
The determination of not only the WF value but also its homogeneity across the wider surface is crucial. Despite the importance of UPS, there has been relatively little discussion on how to conduct measurements on patchy samples with laterally heterogeneous WF. 24−28 At the same time, such heterogeneity is expected in uncoalesced films resulting from island growth, typical of 2D materials. Local WF measurements that can be obtained by Kelvin probe force microscopy (KPFM) are missing for 2D-MoO 3 in ultrahigh vacuum (UHV). In contrast to previous reports, we combine global WF measurements with a nanoscale KPFM study to link the actual functional properties of MoO 3 layers with their atomic structure. Such an approach allows us to point the crystallinity and coalescence of MoO 3 monolayers as important factors for optimization of buffer layers in effective tuning of WF for organic electronic applications.
Our previous work described the possibility of synthesizing 2D nonstoichiometric α-MoO 3 layers directly on HOPG. 15 While studying the morphology, chemical, and electronic structure of α-MoO 3 , we also observed a substantial WF decrease in air-exposed samples by UPS, consistent with the literature reports. 7,8 This work precisely defined the material. However, crucial questions remain open: how the crystallinity of the interface influences the tuning of WF and whether it is an attractive solution for potential applications. This work defines the important and possible general relation between crystallinity of functional materials and optimization of its properties as a buffer layer for energy-level alignment in organic electronics. This study explains what WF the 2D MoO 3 nanosheets have on a local scale when the UHV environment remains unbroken and how to evaluate WF of uncoalesced monolayers with a global technique accurately.
Herein, we report an investigation of the WF of monolayer (1L) MoO 3 on HOPG by a combination of two techniques: KPFM and UPS. We start with presenting high crystalline quality of the grown monolayers and continue with local WF measurements over micrometer-size areas, that is, mapping contact potential difference (CPD) by KPFM in UHV. Next to nanoscale imaging, we provide a commentary on UPS measurement of the WF of MoO 3 monolayers to verify whether this technique is suitable for the determination of WF attributed to 2D islands on a substrate with relatively lower WF. Particular effort is made toward addressing complications associated with double secondary electron cutoff due to insufficient surface coverage. UPS analysis of various 2D-MoO 3 /HOPG lateral area ratios provides information helpful for determining the high WF of one of the materials. Furthermore, density functional theory (DFT) calculations are performed to elucidate the possible underlying mechanism of the influence of oxygen vacancies in MoO 3 monolayers from two directions: top and bottom face.
Driven by an interest in graphene-based electrodes, the presented 2D-MoO 3 /HOPG system is attractive for fundamental study as a model platform to describe the electronic properties of 2D-MoO 3 /graphene heterostructures. Based on the presented complementary results, we prove that 2D-α-MoO 3 is an excellent candidate for replacing the commonly used amorphous functional oxide layers with crystalline ones. We indicate the crystallinity as a highly efficient way of tuning the WF of electrode materials simultaneously with keeping the sub-nanometer thickness limit. The proposed approach allows for uniform modification of electronic properties of devices by a single monolayer of high continuity. Because 2D-MoO 3 can be produced on a centimeter scale (α-MoO 3 ), 16,29−31 our results are potentially easy to scale as needed for industrial applications. In general, our approach carries important implications for characterizing functional interfaces, especially for the purpose of organic electronics. We emphasize that, among global properties and efficiency of the applied materials, the functional layers should be deeply studied on the nanoscale by electronic structure-sensitive techniques, as this is where potential optimization should be looked for. This approach especially paves the way for designing more complex vdW heterostructures. 32,33 ■ RESULTS AND DISCUSSION

Monolayer MoO 3 �Morphological Characterization.
We synthesized monolayer MoO 3 directly on HOPG in UHV through thermal evaporation of MoO 3 powder at elevated substrate temperatures, as described in our previous work. 15 First, we confirm the growth of a monolayer-thin film. Second, their high crystalline quality is verified. These two features play a crucial role in the electronic properties of MoO 3 /HOPG heterostructures, as described in the following sections.
The α-phase of MoO 3 is a thermodynamically stable polymorph with a layered structure consisting of chemically bonded octahedral nets aligned in (010) planes linked by weak vdW forces. 12,22,34 The unit cell of bulk α-MoO 3 contains two layers separated by a vdW gap, with the following lattice parameters: a = 3.96 Å, b = 13.86 Å, and c = 3.70 Å. 17,18,22,34 A monolayer of α-MoO 3 corresponds to half of the unit cell in the b direction, which corresponds to a thickness of ∼6.93 Å, and the surface unit cell is defined by c × a�that is, 3.70 × 3.96 Å 2 . Figure 1a illustrates MoO 3 forming a monolayer (1L) on a graphene-like substrate, that is, HOPG.
We checked the thickness of the grown MoO 3 islands using local conductivity AFM (LC-AFM) under UHV conditions without any exposure to air (see the topography in Figure 1b). The thickness of a layer is presented at the profile (Figure 1b), and an extracted step height of around 6.3 Å confirms a monolayer (1L) growth within experimental error of LC-AFM.
LC-AFM is known to resolve atomic contrast in the current image despite a relatively dull tip and allows for identification of atomic periodicity. 35 The atomic periodicity resolved in the current LC-AFM image in Figure 1c indicates the high crystalline quality of the MoO 3 monolayer. Note that current imaging is possible due to slight nonstoichiometry of the obtained layer, which induces additional states in the band gap. 15 In conjunction with the LC-AFM evaluation of the crystalline quality of the 2D MoO 3 film, high-resolution transmission electron microscopy (HR-TEM) analysis was conducted to assess its crystal structure. TEM-suitable specimens were made by peeling off multilayer graphene flakes with MoO 3 , from the top surface of the bulk HOPG (see the Methods section for more information). The flakes are shown in Figure S1a For the grown α-MoO 3 monolayers, the routinely acquired atomic-resolution scans by LC-AFM together with HR-TEM analysis confirm the high crystalline quality. Next, we focus on the characterization of electronic properties�that is, the WF of the 2D-MoO 3 /HOPG system. Local Surface Potential on Monolayer MoO 3 and HOPG�KPFM Imaging. We employed KPFM to measure CPD, which is directly related to WF of MoO 3 monolayerthick islands on HOPG. For this purpose, the growth process was stopped before obtaining 100% substrate coverage. The asprepared samples possess HOPG regions uncovered by MoO 3 , which serve as a reference during WF measurements. The KPFM measurements were performed in UHV to minimize environmental influence. The UHV conditions are vital because the WF of 2D-MoO 3 is critically influenced by air contaminants: in air, the WF is significantly reduced from 6.9 eV down to ∼5.2 eV for bulk MoO 3 7,8 and to ∼4.7 eV for 2D MoO 3 . 15 The topography image (Figure 2a) shows the MoO 3 monolayer grown on HOPG�a typically shaped island of MoO 3 consisting of multiple domains. 15 Examination of the 1L MoO 3 flake boundaries shows an artificial edge height of around 10 nm; it is an artifact due to uncompensated electrostatic tip−sample interaction and thus appears in regions where the tip apex interacts with both materials having a significant difference in CPD. 36 In this image, one can also recognize that the MoO 3 flake is divided by trenches associated with leaf-like structures on its edges. This division indicates the presence of strong perimeter diffusion toward the tip of growing leaves and some crystallographic misalignment preventing the healing of the trenches. Finally, it is worth pointing out that MoO 3 tends to form a complete monolayer first, and subsequent deposition leads to the formation of islands of multilayers. 15 However, even in uncoalesced monolayer films, some multilayers can be found (∼100 nm of lateral dimension). This observation indicates lower adsorption energy of incoming clusters on edges than on terraces. 15 The corresponding CPD map in Figure 2b shows a clear contrast between the HOPG substrate and MoO 3 . The white color on the HOPG surface is attributed to a high CPD signal, which corresponds to low WF values, while the pink/black color at the MoO 3 monolayer indicates lower CPD and higher WF. The line profile taken over the MoO 3 /HOPG surface in the CPD map (inset in Figure 2c), as indicated in Figure 2b, exhibits two plateaus. The WF difference between the HOPG and 2D-MoO 3 extracted from the profile is negative (ΔCPD MoOd 3 /HOPG = −1.80 V), indicating that the WF of the MoO 3 is higher than that of HOPG.
A histogram extracted from the CPD map is given in Figure  2c. For simplicity, the scale on the right represents WF values as calibrated using WF of HOPG. For quantitative evaluation, data were fitted by two Gaussian distributions, from which mean CPD and full width at half-maximum (FWHM) are determined for both MoO 3 and HOPG: CPD MoOd 3 = −2.10 V (peak #1) and CPD HOPG = −0.30 V (peak #2), where FWHM for both materials is at the same level of 0.03 V. This leads to the surface potential difference between the 1L-MoO 3 islands and the HOPG background of ΔCPD MoOd 3 /HOPG = −1.80 ± 0.05 V. We also reference CPD measurements to a known WF, which allows calculation of WF of 2D-MoO 3 . Thus, with the calibration to the WF value of HOPG, Φ HOPG = 4.50 ± 0.10 eV (see UPS results in the next section), and the MoO 3 2Dislands exhibit WF as high as 6.30 ± 0.12 eV.
In regard to WF homogeneity of the MoO 3 monolayer, it is also worth noting how growth defects affect the WF. Some topographic protrusions can be seen in Figure 2a on the lefthand side. On the corresponding CPD map in Figure 2b, a faint contrast pointing toward slightly lower WF can also be seen besides grain edges, which coincides with the topographic protrusions in the upper left corner. This variation is quantitatively shown in the histogram (Figure 2c) at count columns marked #3 (at −2.00 V) next to the #2 peak. While we assign peak #2 to the area of the smooth MoO 3 monolayer, #3 can be assigned to the lower WF at topographic defects. See also Figure S2 with histograms considering selected areas over the CPD map. Note that the surface potential differences at MoO 3 are within approx. 0.10 V over areas with structural defects.
Furthermore, we observed that the high WF is achieved independent of MoO 3 grain size. For small grains and big grains, a plateau of CPD is observed at the same level�see Figure S3 with multiple grains of various sizes.
Additionally, we support our CPD measurements independent of topographical height variations. Figure S4 shows HOPG steps of height of about 11 nm; the corresponding CPD map displays homogenous CPD distribution, that is, both HOPG terraces are covered by MoO 3 , and are indistinguishable on the CPD maps.
WF of Monolayer MoO 3 /HOPG�UPS Characterization. Next, we use UPS to confirm the local WF with a global method and show how UPS determines WF of the uncoalesced 2D-MoO 3 film. Here, we characterize samples of fractional coverage and the coalesced 1L MoO 3 film. Figure 3 shows three steps in the evolution of the UPS spectrum at the lowkinetic energy region: from pristine HOPG, although partially covered substrates by islands of MoO 3 monolayers, and complete coverage of MoO 3 . The WF is determined from the intersection of linear extrapolation of an edge of secondary electron cutoff (SEC) with the background; 5,15 the x-axis, that is, photoelectron energy, is calibrated to the Au(111) Fermi level (energy�E F ) such that the energy indicated by SEC  We begin by measuring the WF of our reference sample� the UHV-cleaved pristine HOPG substrate. The UPS spectrum in Figure 3a presents an SEC, from which one can read a WF of 4.50 ± 0.10 eV. This value agrees with literature values given in the range of 4.40−4.60 eV. 15,37,38 With fractional coverage by MoO 3 monolayers on top of HOPG (Figure 3b), the SEC significantly shifts to higher energy, indicating a higher WF. However, the high-intensity peak corresponding to the SEC consists of two components observed at 6.15 ± 0.15 and 6.50 ± 0.15 eV. We assign the first SEC at low energy (highlighted purple) to an apparent WF value of both: uncovered HOPG patches (blue) and top MoO 3 islands (red), while the second SEC at high energy (highlighted red) corresponds to the MoO 3 monolayer only. The double SEC can be explained by electrostatic potential distribution above a two-component surface. The electrostatic potential of MoO 3 , that is, the high WF patches, creates an additional energy barrier for photoelectrons emitted from uncovered HOPG, that is, the low WF patches. The explanation of the observed effect is well described by Schultz et al. 25 More on this topic will be covered at the end of this section. Details on fitting distinct SECs as step edges in CasaXPS are presented in Figure S5.
Finally, after depositing a coalesced MoO 3 monolayer film (Figure 3c), the SEC interpretation is again straightforward, and the WF of MoO 3 approaches 6.40 ± 0.10 eV in agreement with our KPFM measurements. This confirms that the highenergy SEC in Figure 3b provides an accurate measure of 2D MoO 3 WF (see also Figure S6 with overlaid spectra for direct comparison). The WF values we obtained are summarized in Table 1 together with literature values for comparison.
It should be noted that the shape of the SEC region is a measure of both uniformity of WF and surface cleanliness. 27,40 Again, the width of SEC increases with inhomogeneity of the WF within the probed area: the width of double SEC of fractional coverage by 2D-MoO 3 (Figure 3b) is wider than single SEC of the coalesced film (Figure 3c). The monolayer of MoO 3 features homogenous WF (the SEC width, i.e., FWHM = 0.05 eV), as shown in Figure 3c. Moreover, since air-derived contaminants lower the WF of clean MoO 3 , 15 the absence of any low kinetic energy tail indicates high surface cleanliness (Figure 3c).
At this point, it is worth discussing the interpretation of UPS for studies of uncoalesced films or 2D flakes with higher WF than their substrate. UPS has become, in recent years, a popular technique for measuring the efficacy of WF tailoring for photovoltaic and light-emitting devices. [2][3][4][5]8,24 Our setup seems to be a perfect system to discuss the limitation of this technique for nonhomogenous systems. As mentioned earlier, in the literature, some efforts focus on patchy samples with laterally nonhomogeneous WF, where patchiness derives from dirt or regular patterning. 24−28 Here, we present a real-life example where area of <1 nm-thick patches is irregular and depends on the morphology of the substrate because MoO 3 tends to nucleate at terrace edges of HOPG. The topography of our intentionally prepared test sample can be described as randomly distributed MoO 3 islands, which take different shapes, but it is worth noting that there is a gradient of surface coverage from dense on one side to scarcer on the other (Figure 4a). For visualization of the coverage gradient, additional photography was taken on an interference pattern of the SiO 2 substrate with the MoO 3 thin film�see the inset in Figure 4b. Due to this gradient, we can investigate different coverages of MoO 3 over the underlying HOPG within the probed area by varying the position for the UPS measurement across a sample with the uncoalesced monolayer film. We took a range of spectra at various points (see lines in Figure 4a) to show how acquired UPS spectra can vary under such conditions. As expected, the UPS spectra are consistent within the horizontal line, where coverage is similar�see Figure 4b. On the other hand, when the MoO 3 /HOPG area ratio varies�vertical line in Figure 4a�the UPS spectra significantly change. Figure 4c presents the SEC region of that series of spectra, while Figure 4d replots them with normalized intensity to highlight the positions of high-energy SECs, corresponding to the WF value of MoO 3 , which are marked by dotted vertical lines. The shift in the low-energy SECs is clear (about up to 200 meV), and the change in SEC shape, width, and relative intensities of distinct SECs is believed to depend on the detail of the MoO 3 /HOPG morphological arrangement within the probed area. Regarding the latter, note that if the area covered by MoO 3 is insufficient to contribute and build up a distinct high-energy SEC feature, a spectrum reveals only one SEC feature�see the top blue spectrum in Figure 4c. In this case, the pure MoO 3 contribution is hidden, and the single broad SEC reveals a WF higher than in bare HOPG and lower than that of MoO 3 due to a higher electrostatic barrier for photoelectrons coming from HOPG as a result of the presence of MoO 3 islands. To discuss the usability of UPS for WF determination of 2D islands, some useful information can be deduced from these measurements, as listed below: (I) for different coverages of the sample (in the form of 2D-MoO 3 islands), the high-energy SEC (if visible) is uniformly 25 located in our case at ∼6.50 eV; thus, it is determined as the absolute value of WF of a higher WF material, that is, MoO 3 ; (II) the energy of the low-energy SEC shifts; thus, it is dependent on the area ratio 25 between the substrate and high WF material (here, exposed HOPG and MoO 3 ) because the position of low-energy SEC is governed by electrostatic potential distribution above a two-component surface; (III) with the low-energy SEC shifting toward lower energies, the relative heights of the peaks change, that is, the contribution from HOPG increases, and thus, the intensity of the low-energy SEC increases; likewise, the overall intensity of a spectrum in the SEC region also increases 25 �the same variation in SEC was found for patches composed of ITO and Au by Sharma et al.; 26 and (IV) to observe a distinct high-energy SEC, a minimum area of a high WF material is required. Below this coverage threshold, the UPS measurements do not allow for clear conclusions about effective electronic properties of a material. DFT Calculations of Electrostatic Potential of Monolayer MoO 3 . The experimental finding of the WF of pristine monolayer α-MoO 3 grown on HOPG approaches 6.40 eV, which is lower than bulk α-MoO 3 (6.9 eV 2,6−8 ). The theory predicts that the electronic structures of a monolayer and bulk α-MoO 3 do not alter significantly. 12,16−19 Nevertheless, the actual WF is sensitive to many factors, including the underlying substrate and stoichiometry. 6,36,41 The latter is known to influence the WF value of TMOs. 42−44 Moreover, the possible LC-AFM current image shown in Figure 1b may indicate a slight nonstoichiometry in our films. Even a slightly nonstoichiometric α-MoO 3 monolayer is conductive; thus, it may not reach the highest possible WF for the bulk material. Here, we present theoretical predictions of WFs of stoichiometric and nonstoichiometric monolayers.
The WF can be estimated by the DFT calculations, based on electrostatic potential across the interface of a material and vacuum. The calculations were carried out for two monolayer structures: fully oxidized MoO 3 and reduced MoO 2.875 (with a single oxygen vacancy per 2 × 2 supercell, located in the bottom layer; see Figure 5).
The electrostatic potential predicted by our DFT calculations averaged over the monolayer (xy) plane is plotted in Figure 5 as a function of the z coordinate which is orthogonal to the monolayer plane (along the b axis), that is, across the monolayer−vacuum interface. The zero value of the z coordinate is selected in the plane in between Mo atoms (precisely at the average of z coordinates of all the atoms for the studied structure). The location of the atoms along the z axis in the calculated unit cells is marked in the plot with vertical lines (see Figure 5). For the MoO 3 monolayer without defect (marked red), the distance dependence of the electrostatic potential is fully symmetric with respect to the middle plane. Two deep potential minima correspond to Mo atoms and pairs of nearby O atoms, whereas much shallower minima mark the positions of outer O atoms in the structure. The WF value amounts to 6.67 eV.
The symmetry mentioned above is perturbed for the MoO 2.875 monolayer (see the blue plot in Figure 5) with a single O vacancy. For the upper side with no vacancy (z > 0), the increase in the interatomic distances along the z axis is the only noticeable feature (see the shifted minimal positions and atomic coordinates). The predicted WF value on the defectfree side is only slightly increased with respect to the stoichiometric MoO 3 and amounts to 6.72 eV. At the bottom side with the vacancy defect (z < 0), the atomic positions are much more perturbed, as visible in shallow and manifestly asymmetric potential minimum, also shifted due to increasing interatomic distances. The WF predicted at the defect side is reduced to 6.36 eV. Such reduction suggests that our experimental data are affected by imperfections in the stoichiometry of MoO 3 . Our 2D-MoO 3 was likely slightly oxygen-deficient, in line with the previous reports on thermally evaporated MoO 3−x . 8, 15,45 So far, we have focused on the WF estimated from the constant values of the planar average electrostatic potentials at a considerable distance from the surfaces. We additionally refer the reader to Figure S7, which shows WF maps for both structures in certain z distances over the surfaces�it is particularly interesting how WF is distributed over the MoO 2.875 monolayer from each side. These WF maps theoretically predict the ultimate resolution KPFM measurement, which cannot be achieved in real-life experiments.

■ CONCLUSIONS
The MoO 3 monolayer films presented here exhibit electronic properties, that is, high WF, similar to bulk MoO 3 , while being <1 nm thick and quasicontinuous. We achieved a WF of approximately 6.40 eV for the crystalline MoO 3 (010) monolayer on HOPG(0001); this value was reproducible across the samples, provided that UHV conditions were used for both their growth and characterization. The DFT calculations predict the significant influence of nonstoichiometry on the WF value and indicate that the oxygen vacancies on one side of the MoO 3 monolayer cause the variation in WF of an opposite sign for each face.
In contrast to the existing reports, we emphasize that crystalline 2D-MoO 3 guarantees high WF, and in this regard, the material does not require thicker layers but coalesced monolayers. With maturing synthesis of continuous monolayer films, which tend to wet graphene-like substrates, 2D MoO 3 becomes a great candidate for a monolayer buffer layer for WF enhancement.
Our observations not only address the critical issue of the WF for 2D MoO 3 under UHV conditions but also represent a real-life example of how the commonly accessible UPS apparatus can be used to determine WF attributed to 2D islands or flakes on a substrate with lower WF. Thus, our study showcases a simpler way to obtain a WF of 2D nanomaterials using UPS only, without the need for additional characterization techniques, which may also act as a tool supporting future investigations. These observations carry important implications for the characterization of 2D materials beyond MoO 3 and the design of vdW heterostructures for energy-level alignment with graphene-like substrates. ■ METHODS Monolayer MoO 3 Synthesis. The synthesis was the same as that used in ref 15, namely, films of α-MoO 3 were thermally evaporated on heated to 220°C HOPG(0001) substrates (ZYA grade) from the powder source of MoO 3 (Sigma-Aldrich). Freshly deposited films were loaded into a vacuum suitcase and transferred to the analysis chamber of an ultrahigh vacuum (UHV) Multiprobe P (Scienta-Omicron) system with a base pressure of 5 × 10 −10 mbar, where all UPS and AFM measurements were performed at room temperature.
UPS Characterization. A hemispherical energy analyzer Phoibos 150 (SPECS), with a 2D-CCD detector, was used for photoelectron spectroscopy measurements. UPS studies were performed with the helium I line (21.22 eV) of VUV Source HIS 13 (Focus GmbH) on samples perpendicular to the analyzer and with an applied bias of −3.1 V. 27 At a pass energy of 5 eV, the analyzer offers an energy resolution of ∼130 meV as determined from the Fermi edge of Au(111) for the settings used in this study (step size of 0.03 eV). The calibration was conducted at the Fermi level of Au(111).
Analysis of the photoemission spectra was performed using CasaXPS software. The WFs were derived from the intersection of linear extrapolation of a steep edge of SEC with the background, that is, edge-up background type provided in CasaXPS (see Figure S5). AFM Characterization. Local-conductivity AFM (LC-AFM) and frequency-modulated KPFM (FM-KPFM) experiments were performed with an Omicron Matrix system and a Nanonis controller (SPECS) with bias applied to the tip during the measurement. CPD maps were obtained simultaneously with topography noncontact-AFM (nc-AFM) images. The probes were highly doped silicon NANOSENSORS PPP-NCH with a nominal resonance frequency and force constant of 330 kHz and 42 N/m, respectively, with an AFM tip radius of curvature of <10 nm when first used. The scanning conditions were chosen to highly reduce the tip−sample distance while maintaining the stability of measurement, and the frequency shift setpoint near −20 Hz was used. This allows us to improve spatial resolution of CPD and minimize topographic artifacts at MoO 3 edges, which were a significant issue due to the very large WF difference between the materials.
The contrast of a CPD map is directly related to WF differences on the surface. The WF of MoO 3 was calculated with respect to the WF of HOPG measured using UPS (HOPG was cleaved in UHV before the measurement; see Figure 3a) according to the following equation: Φ MoOd 3 = Φ HOPG + e × (CPD HOPG − CPD MoOd 3 ), where e is the elementary charge. The KPFM images were processed using Gwyddion software. 46 TEM Measurements. TEM imaging was performed at the National Center for High Resolution Electron Microscopy (nCHREM) at Lund University. The MoO 3 /HOPG samples were a few weeks old and subjected to exposure to an inert (Ar) atmosphere for transportation and laboratory air (typically below 5 min) to prepare TEM-suitable specimens. The MoO 3 /HOPG flakes were peeled off and mounted on a TEM grid (lacey formvar-carbon film on 200 copper mesh, Ted Pella, Redding, USA). The imaging was performed in a JEOL JEM-2200FS TEM at 200 keV in a low-dose mode (total dose < 200 e Å −2 ), thus minimizing possible beam damage. The images were recorded on an F416.0 camera (TVIPS), using Serial EM software. 47 Computational Methods. The ab-initio calculations were performed using the QUANTUM ESPRESSO suite, implementing DFT 48 using the plane-wave pseudopotential method. 49,50 The scalar relativistic pseudopotential approach with the projector augmented wave method 51 was used, with the Perdew−Burke−Ernzerhof exchange correlation functional 52 accepted. The empirical van der Waals correction 53 and the dipole correction 54 were taken into account. The 2D layers were modeled in slab geometry; the calculations involved a monolayer of MoO 3 and a monolayer of MoO 2.875 modeled using a 2 × 2 supercell with a single O atom removed from the lower atomic plane. In all the calculations, the total cell size in the direction perpendicular to the layer was set to 50 Å to ensure sufficient vacuum. For structural optimization, the variation of both in-plane lattice constants and all the atomic positions was allowed, ruled by the Broyden−Fletcher−Goldfarb−Shanno (BFGS) quasi-Newton algorithm, based on the trust radius procedure. All the structural relaxation and self-consistent calculations were based on 24 × 24 mesh of k-points. For the monolayer MoO 3 structure, the convergence criteria were as follows: 10 −11 Ry for energy and 10 −9 Ry/Bohr for force. For MoO 2.875 , the convergence criteria were as follows: 10 −6 Ry for energy and 10 −4 Ry/Bohr for force.

* sı Supporting Information
Additional HR-TEM images; KPFM images with detailed histogram analysis; UPS spectra of MoO 3 / HOPG including fitting of step edges at SECs and direct comparison of the discussed spectra; and DFT calculations of electrostatic potential maps of MoO 3 and MoO 2.875 monolayers (PDF) ■ AUTHOR INFORMATION Corresponding Authors